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Abstract 

We define a notion of logarithmic, Coulomb and Riesz interactions in any dimension 
for random systems of infinite charged point configurations with a uniform background of 
opposite sign. We connect this interaction energy with the “renormalized energy” studied 
by Serfaty et al., which appears in the free energy functional governing the microscopic 
behavior of logarithmic, Coulomb and Riesz gases. Minimizers of this functional include 
the Sine-beta processes in the one-dimensional Log-gas case. Using our explicit expression 
(inspired by the work of Borodin-Serfaty) we prove their convergence to the Poisson process 
in the high-temperature limit as well as a crystallization result in the low-temperature limit 
for one-dimensional systems. 


1 Introduction 


1.1 General setting 

Logarithmic, Coulomb and Riesz interactions. We consider a system of points (which 
we can think of as being point particles carrying a positive unit charge) in the Euclidean space 
interacting via logarithmic, Coulomb or Riesz pairwise interactions 


(1.1) 

g{x) 

= — log x|. 

in dimension d = 1 

(1.2) 

9{x) 

= -logkl, 

in dimension d = 2 


or in general dimension 


(1.3) 



max(0, d — 2) < s < d. 


Cases dni) and (na) are known as one- and two-dimensional log-gases, and we will refer 
to them as the “logarithmic cases”. One-dimensional log-gases have been extensively stud¬ 
ied for their connection with important random matrix models known as the /3-ensembles 
(see |ForlO] i. The two-dimensional log-gas is known in the physics literature as a two-dimensional 
one-component plasma (see |A.T81j l and can also model non-Hermitian random matrices such as 
the Ginibre ensemble |Gin65j . The cases (11.31) correspond to higher-dimensional Coulomb gases 
(if s = d — 2) or Riesz gases. 

The statistical mechanics of N points (xi,..., xn) interacting pairwise via g under a confining 
potential V at inverse temperature /3 G (0, -|-oo) is given by the canonical Gibbs measure 


dFjv/3(xi,... ,xjv) := ,,, ^xn, 
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where is a normalizing constant. The macroscopic behavior i.e. the behavior of the 

empirical measure /tat := ^xi is well studied in the limit N —)• oo, see e.g. |Serl4j and 

the references therein. The limiting macroscopic arrangement is described by the equilibrium 
measure Heq, which is a probability measure on depending on V, with compact support S, 
such that {/tatIat converges weakly to peg, Tw./S-a.s. 

In order to study the microscopic arrangement of the particles, Sandier-Serfaty have derived 
in |SS12j (see also |SS14j . [E,S15j . |PS15j l a second-order energy functional Wn which governs 
the fluctuations around peg, together with an object W defined on infinite point configurations 
which is the limit of VFat as —)• oo in the sense of T-convergence of functionals (see [Serl4| l. Let 
x' := u’jq := and let p'^^ be the push-forward of peq by x The “blown- 

up” point configuration encodes the position of particles at the microscopic (inter-particle) 
scale We let 

(1.4) ITAr(xi,... ,XAr) := ^ ~ “ dpiqix)) 0 (du'^ - dp'^g{y)), 

where A is the diagonal. It can be seen that ILAr computes (up to the factor ^) the Coulomb/Riesz 
interaction of the electric system made of the finite (charged) point configuration u'jq and of a 
negatively charged background of density dp'^^, with itself, without the infinite self-interactions 
of the charges because the diagonal A is excluded. 

An integration by parts shows that Wm may be re-written with the help of the associated 
electric field :=Vg* {v'jq — p'(,q), whose norm is computed in a renormalized fashion to take 
care of the singularities around each charge (we will come back to this procedure in Section 13.ip . 

Following the same procedure, a renormalized energy functional W is defined on the space of 
electric fields corresponding to infinite point configurations together with a uniform background 
of intensity 1. If C is a point configuration, its Coulomb/Riesz energy W(C) is then defined as 

W(C) := inf>V(F;), 

E 

where the infimum is taken among the set of electric helds E which are compatible with C. 
Finally if P is a random point process (a probability measure on point configurations) its energy 
is defined by 

Wclec(p) _ Ep[W]. 

We refer to Section [3.II for more details. The superscript “elec” is added by us and refers to this 
“electric” approach to the definition of a Coulomb/Riesz energy. 

Free energy at microscopic scale. In [LS15j S. Serfaty and the author have obtained a 
second-order (or process level) large deviation principle concerning the average microscopic be¬ 
haviour of the particles under the canonical Gibbs measure at inverse temperature /3 € (0, -|-oo). 
This behaviour is characterized by a certain random point process P (the law of a random point 
configuration) and it amounts to minimizing a free energy functional of the form 

(1.5) Pf 3 {P) := /3W®'®‘=(P) + ent[P|n] 

on the space 'Psq(T) of translation-invariant random point processes whose mean density of 
points is 1. The term ent[P|n] denotes the specific relative entropy of P with respect to the 
Poisson point process of intensity 1 on it is the infinite-volume analogue of the usual relative 
entropy. 

From a physics perspective, knowing the minimizers of the free energy and how they behave 
as (3 varies allows one to retrieve some of the thermodynamic properties of the physical system at 
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the microscopic scale (e.g. the existence of phase transitions). From the random matrix theory 
point of view, it was proven in |LS15j that an important family of point processes governing the 
microscopic behavior of eigenvalues, namely the Sine^j processes of Valko-Virag (see |VVn9j i. 
minimizes for /3 > 0. Hence it would be very useful to have information on J-p, its level sets 
and its minimizers (depending on /3). A drawback of the free energy Tp is that computing it 
explicitly is hard. The energy term in particular is difficult to evaluate and except for the case 
of periodic configurations (for which exact formulas hold) no value of W is known and the mere 
finiteness of W(C) for a given point configuration C is unclear in general (see however |(fS13j for 
some criteria). On the other hand, level sets of W are easily seen to be degenerate because small 
perturbations of any given configuration C typically do not change its energy. Concerning both 
issues, it turns out to be helpful to look for a definition of the energy directly at the level of 
stationary random point processes instead of averaging the energy computed configuration-wise. 


The approach of Borodin-Serfaty. In |BS13j a related notion of a renormalized energy 
for random point processes was introduced in the logarithmic setting (jl.ip and Given 

a stationary random point process P, Borodin and Serfaty proceed by periodizing the point 
process induced in a square (or interval) of sidelength R and computing its renormalized energy 
by the mean of explicit formulas valid in the periodic setting. If (ai,... ,aAr) is a configuration 
in a torus T of volume N in the associated periodic point configuration has an energy 


W(ai,... ,aAr) 





ttj) 


+ cig lim 
x-s-O 


(gp'^''(x) + 


log^\ 

Cd,s J 


where Cd,s is a constant and is a certain periodic Green function which has a logarithmic 
singularity at 0. Taking the expectation under P, using an expansion of GP®’’ and sending 
N ^ oo they obtain an energy W®®(P) (our notation) which may be written, up to an additive 
constant, as 

( 1 . 6 ) W^^{P)= f -log\v\{p2,p{v)-l)dv, 

where p2,p denotes the two-point correlation function of P, which can be seen as a function of 
one variable by stationarity (we abuse notation and let P2,p{v) '■= P 2 ,p{ 0 , v)). Using this explicit 
expression in terms of P 2 ,p, they are able to compute the energy for some specific point 
processes (e.g. the Sine^ processes for /3 = 1,2,4, and the Ginibre point process) as well as 
to solve minimization problems (the minimization of over a large class of determinantal 
point processes). However no general rigorous connection is drawn between W®® and the electric 
definition W®^®® which derives from the energy functional Wat. Moreover the formulas of |BS13] 
only apply to random point processes for which P2,p{x, y) — l decays fast enough as \x — y\ —)• oo. 
The approach of the present paper is strongly inspired by the one of [BS13] and is an attempt 
to give a partial connection between and W®*®®. 


1.2 Main results 

The purpose of this paper is twofold. First we introduce an energy (“int” as “intrinsic”) 
defined on Vsp{X) which is expressed only in terms of g and of the two-point correlation function, 
and we connect W™** with W®^®®. In a second part we use W™* to handle the energy term in Pp, 
which allows us to describe the behavior of minimizers of Pp in the limiting cases /3 ^ 0 (in any 
dimension) and /3 —)■ oo (in dimension I). 
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New definition of the Coulomb/Riesz energy. If P is a stationary random point process 
of intensity 1 and p 2 ,p denotes its two-point correlation function, we define in Section 13.21 its 
“intrinsic” energy (with respect to the interaction g) as 


W“t(P) —liminf 


// g{x,y) {P2 {x,y)- l)dxdy, 

JJcl\A 


where Cp is the hypercube [—■§, An equivalent formulation is 


(1.7) W‘"*(P) := liminf ^ f 9iv){p2,p{v) - 1) ]J(P - \vi\)dv, 

R^oo J[-R,R]d\{0} 


where v = (ui,... ,Vd) and where we made again the abuse of notation P2,p{v) ■= P 2 ,p{ 0 ,v). 
The expression (|1.7I) shows similarities with (ll.bp in the logarithmic cases. 

Let us also define, in the logarithmic cases 


(1.8) P'°§(P) := lim sup i^ff^ iP2,pix, y) “ ^)dxdy -h 1 j log P, 

where (7^°® is a constant whose value is irrelevant for our concerns. 

Finally we introduce the free energy functional analogous to Tp (defined in (ll.5|) l 

;= /3W“* + ent[-|n], 

or in the logarithmic cases 

T'p := /3(W“* + + ent[-|n]. 

Let us recall the following definition: let A be a topological space and f,g : X —>• M two 
functions. We say that g is the lower semi-continuous regularization of / if for any x € X we 
have 

g{x) = lim inf/(y). 

y^x 

Our first main result is 

Theorem 1. The functionals and are related as follows. 

• In the one-dimensional logarithmic case (ll.ip . is the lower semi-continuous regu¬ 
larization o/-I- and for any (3 G (0,-|-oo), Xp is the lower semi-continuous 

regularization of X'^. 

• In the non-Coulomb cases (m with s > d — 2, is the lower semi-continuous regular¬ 

ization o/W™*, and for any j3 G (0,-|-oo), Xp is the lower semi-continuous regularization 
ofX'^. 

• In the two-dimensional logarithmic (Coulomb) case ()1.2p . we have 

Welec < 

• In the higher dimensional Coulomb cases (m with s = d — 2, we have 

W®'®® < W“^ 
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A first interest of Theorem [T] is that it provides a way of showing that a given random point 
process has finite energy. For example in the d = 3 Coulomb case, the Poisson point process of 
intensity 1 is easily seen to satisfy W^“*(II) = 0, hence is finite and nonpositive. 

Let us emphasize that Theorem [1] is less precise in the Coulomb cases than in the cases (|l.ll) 
and do]) with s > d — 2, to which we will henceforth refer as the “non-Coulomb cases”. 

In the following statement, by saying that two minimization problems are equivalent we 
mean that both functionals have exactly the same infima. If g is the lower semi-continuous 
regularization of / on X, then the minimization problems associated to / and g are equivalent, 
thus 

Corollary 1.1. We deduce from Theorem{l\ that in the non-Coulomb cases 

1. Minimizing onVs,i{W) is equivalent to minimizing (orW^^^onVs,i{X). 

2. For any (5 € (0,-|-oo), minimizing Tg on Vs,i{W) is equivalent to minimizing on 
Vs,i{X). 

Applications. We use the explicit expression of W™* and its link with to perform simple 
computations, which allow us to prove the following two results concerning the minimization of 
Tg in the limit /3 —)• 0 and /I —oo. 

High-temperature limit. As can be expected, when /3 —)• 0 any minimizer of J-g gets 
close to the minimizer of the entropy term. 

Theorem 2. For all cases dm), dn, dn, the minimizers of converge as 13 ^ 0 to the 
law of the Poisson point process H. Moreover this convergence holds in entropy sense i.e. 

(1.9) lim sup ent[P^|n] = 0. 

In the special case of one-dimensional log-gases, as proven in |LS15j a minimizer of Tg is the 
Sine /3 process of Valko-Virag [WOhj . Hence our method yields another proof for a recent result 
of Allez and Dumaz |AD14j : 

Corollary 1.2. As /3 —>• 0 the Sineg point process converges weakly in the space of Radon 
measure (endowed with the topology of vague convergence) to the law of a Poisson point process 
on M. 


Low-temperature limit. In dimension 1 we may also characterize the limit ft ^ oo (the 
low temperature limit) of the minimizers of Pg. We let Pz be the stationary random point 
process associated to the lattice Z 


( 1 . 10 ) 


•— / dx-\-z dx, 


which can also be seen as the law of the point configuration u-\-Z where u is a uniform random 
variable in [0,1] and where we let x -|- Z denote the point configuration {x k,k € Z}. 


Theorem 3 (Crystallization for d = 1). For d = 1 and in both cases (II.ip or (jl.hp . the random 
point process Pz is the unique minimizer o/W®^®® on Vsp{W)- Moreover if {Pg}g is a family of 
minimizers of we have 


( 1 . 11 ) 


lim Ps = 
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Theorem [3] is a crystallization result, proving the convergence to the one-dimensional crystal 
as /? —)• oo. A similar result was proven in the one-dimensional logarithmic case in |Lebl5] by 
using the explicit expression available in the periodic setting together with an approximation 
by periodic point processes. The method here is similar in spirit but follows a simpler approach 
which works for Riesz cases as well. 

1.3 Open questions 

Let us now briefly mention some questions that are raised by, or related to our present study. 

Minimizers of W™*. Could one determine the infimum of W™* (thus of in the non- 

Coulomb cases) thanks to its explicit form? A necessary condition on the realizability of p 2 as a 
two-point correlation function of some stationary random point process of intensity one is that 
T 2 ■= P 2 — ^ satishes 

T 2 > -1, Tb > -1, 

(see |KLS07[ Section 2]) where T 2 denotes the Fourier transform of T 2 (in a sense that should, 
in general, be precised). Thus we may start by asking whether the linear optimization problem 
of minimizing 

liminf [ T 2 {v)g{v)Y[ dv 

R -^00 R J 

can be solved on the convex set {T 2 > —1,T2 > —1}. It is unclear to us whether the symmetry 
of the constraints on T 2 and T 2 might be of any use (however let us observe that the two-point 
correlation function of Pz is its own Fourier transform, and that the expected minimizer in 
d = 2, namely the stationary random point process associated to the triangular lattice, exhibits 
a self-duality of the same kind). 

Decorrelating random point processes. We may also investigate the problem of minimiz¬ 
ing over particular sub-classes of We have already mentioned the work of |BS13j 

where this minimization is considered over determinantal point processes in dimension d = 1,2. 
Another interesting aspect is that of “decorrelating” random point processes. Let us say that P 
is a decorrelating random point process (or that P decorrelates) when T 2 {v) := P 2 {v) — 1 tends 
to 0 as |u| —>■ 00 , with some speed criterion to be fixed, or when T 2 lies in some reasonable 
class (e.g. L^-spaces). It is unclear to us whether there is a minimizing sequence for W™* made 
of random point processes which decorrelate. A negative answer would hint at a soft kind of 
phase transition as (5 varies. Indeed as /3 ^ 0 the minimizers of Pp converge to the law of 
a Poisson point process 11 (according to Theorem [2]) , which is a typical decorrelating random 
point process, and as /3 —>■ 00 they would rather leave this class. 

The Id Log-gas case. Such a transition may be formally observed in the one-dimensional 
logarithmic case. Let us recall that for d = 1, s = 0 there exists a “concrete” family of minimizers 
for Pp stemming from Random Matrix Theory (RMT), namely the Sine^ processes. They arise 
as the N ^ 00 limit of microscopic point processes observed in the so-called “Gaussian (3- 
ensembles”, which are a generalization for any (3 of the classical Gaussian ensembles of RMT 
(see |DEn2j and jVV09]). It is known (see [Na kl4] i that Sine /3 is also the law of the limiting 
microscopic point process for the “Circular /3-ensemble”, which is another RMT model. The 
N ^ 00 limit for two-point correlation function of the Circular /3-ensemble has been derived in 
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|For93j for even values of f5. Equation |For93[ (11a)] should thus describe the large-x asymptotics 
of the two-point correlation function of Sine^, i.e. formally we have 

/ 9/2 

T2,Smepix) := P2,Sinep{x) - 1 ~ fkix)x~^^^ ^^, 

k=l 

where /fc(x) is a slowly oscillating function. The leading term of T 2 ^sinep{x) as x —)• oo is of order 
we would thus expect T 2 ^sine^ to leave the class LP as soon as /3 > 4p. 

Decorrelating random point process of minimal energy. A negative answer would 
also raise the question of finding the minimizer (or a minimizing sequence) for W™* among 
decorrelating random point processes. It seems to us that a good candidate for a lower bound 
on the energy is given by the hypothetical “hardcore Poisson point process” Ilhc whose two-point 
correlation function would be P 2 ,hc = 1 ~ Is where B is the ball of center 0 and unit volume in 
In dimension d = 1 it is not hard to construct a sequence of random point processes whose 
energies converge to the associated energy log juj. In arbitrary dimension, it easy to see that 
any sub-Poissonian random point process (i.e. a random point process such that p 2 < 1, e.g. 
any determinantal point process) has a larger energy than Ilhc. 

Plan of the paper and ideas of proof. In Section [2] we give some general dehnitions and 
notation. 

In Section [3] we recall the dehnition of the renormalized energy in the sense of |SSI2j . |RSI5j . 
[PSI5] . then we introduce the alternative object W™* dehned on the space of random point 
processes, and in the stationary case we give a simple expression of W“*(P) in terms of the 
two-point correlation function of P. 

In Section 0] we give some preliminary results. In particular we observe that while is 

by dehnition computed in terms of the energy of global electric helds dehned on the whole space 
]^d+fc, object W™** is rather a limit as ii —)• oo of the energy of local electric helds dehned on 
hyper cubes of sidelength R. 

Section [5] is devoted to the proof of Theorem [TJ The proof goes in two step: hrst we show 
that (or in the logarithmic cases) on the space 'Ps^i(T), then conversely 

(in the non-Coulomb cases) for any P £ Ps^i{X) we prove the existence of a “recovery sequence” 
{Pn}n converging to P and such that limAr^oo W™*^(PAr) < W®^®®(P). In both steps the key 
element is the screening lemma of [PS 1 5] (following |SSI2| . |SSI4] . |RSI5p which, heuristically 
speaking, allows us here to construct a global electric held from a local one, and vice versa. 

In Section 0] we prove Theorem 0] about the convergence to the law of the Poisson point 
process 11 of minimizers of as /3 —>• 0. If W®^®®(n) is hnite we may see directly that ent[P,g|n] 
must go to zero as ,8 —>■ 0, and then the specific Pinsker inequality implies that —)• 11 as 

/? —)• 0. However in some cases the hniteness of W®^®®(n) is false (e.g. d = I, s = 0, see |LSI5p 
or yet unknown (d = 2, s = 0). We use the fact that W®*®® < W™** to construct a sequence 

of random point processes converging to 11 in entropy sense and whose renormalized energy is 
hnite. 

In Section [7] we restrict ourselves to the one-dimensional cases dni or and we prove 
the crystallization result of Theorem [3l We start by using a convexity argument to show that 
Pz is the only minimizer of W™* over More precisely we obtain a quantitative bound 

below on W™*(P) — W“*(Pz) in terms of the two-point correlation function of P. This translates 
into a bound below for W®^®®(P) — W®^®®(Pz) which implies that Pz is also the only minimizer of 
W®^®®. Moreover as ^ oo we show that W®^®®(P,g) must go to W®^®®(Pz), which in turn implies 
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that as /? —)• oo the two-point correlation function of P /3 converges to P 2 ,P^ in the distributional 
sense. Thanks to the “rigidity” of the lattice it is not hard to deduce that in fact P /3 converges 
to Pz as j3 ^ 00 . 

Aknowledgements. The author would like to thank his PhD supervisor, Sylvia Serfaty, for 
helpful discussions and many useful comments on this paper. 

2 Definitions and notation 

2.1 Generalities 

If {X,dx) is a metric space we let V{X) be the space of Borel probability measures on X and 
we endow V{X) with the distance 

(2.1) dp(^)(Pi,P 2 )=sup|y’P(dPi-dP 2 )| PGLipi(X)|, 

where Lip^(A) denotes the set of functions P : A —)• M that are 1-Lipschitz with respect to dx 
and such that ||P||cxd < 1- It is well-known that the distance d-p(^x) metrizes the topology of weak 
convergence on V{X). If P S P{X) is a probability measure and / : A —)• a measurable 

function, we denote by Ep [/] the expectation of / under P. 

The ambient Euclidean space of dimension d is denoted by We will often need to work 
in ^d+k /j = 0 (in the case of Coulomb interactions) or 1 (see Section 12 .4p . 

For any P > 0 we denote by Cr the hypercube [—P/2, P/2]'^ C and by Cr the hypercube 
[—P/2, P/2]'^+^ C Similarly we let Br be the ball of center 0 and radius P in and Br 

be the ball of center 0 and radius P in 

If A is a set, we let A := {(x, x), x G A} C A x A be the diagonal of A. 

2.2 Point configurations and random point processes 

Point configurations. If d is a Borel set of we denote by X(A) the set of locally finite 
point configurations in A or equivalently the set of non-negative, purely atomic Radon measures 
on A giving an integer mass to singletons (see |DV.I 88 j l. We will often write C for 
We endow the set T := T’(M'^) (and the sets X{A) for A Borel) with the topology induced 
by the topology of weak convergence of Radon measure (also known as vague convergence or 
convergence against compactly supported continuous functions), these topologies are metrizable 
and we fix an arbitrary compatible distance. 

The additive group acts on X by translations if C = {xi,i G /} G T we let 

( 2 . 2 ) OfC 

We denote by Mr : T i—>• N the number of points of a configuration in the hypercube Cr, 
and by Vr the discrepancy Vr = Mr — R^. 

Random point processes. A random point process is a probability measure on X. We 
denote by Vs{X) the set of translation-invariant (or stationary) random point processes. We 
endow V{X) with the topology of weak convergence of probability measures on A. A compatible 
distance on Vs{X) is defined in ( 12 . 11 ) . 




Remark 2.1. Another natural topology on V{X) is the “convergence of the finite distributions” 
IDVJOSi. Section 11-1], also called the “convergence with respect to vague topology for the counting 
measure of the point process”. The two notions of convergence coincide as stated in jPVJOtA 
Theorem 11.1.VII]. 

We will use several times the operation of averaging a random point process over translations 
in some measurable subset. If P is a random point process in and K a measurable subset 
of with finite, non-zero Lebesgue measure, we define the average Pffi of P over translations 
in K as the law of the random variable uk + C where uk is uniformly distributed according to 
the normalized Lebesgue measure on K, and C has law P. The sum of a vector and a point 
configuration is defined as the point configuration 

X + C := {x + p,p G C}. 

Hyperuniformity. Following [TSOSj we say that a random point process P is hyperuniform 
if we have 

(2.3) Ep[P|] = 0{R'^-^). 

For a stationary one-dimensional random point process P, hyperuniformity is easily seen to 
be equivalent to the following property: for some r > 0, there are P-a.s. between k — r and k + r 
points in any interval of length k. 

2.3 Correlation functions 

Let P G V{X) be a random point process. For any re > 1 the re-point correlation function 
Pn,p is the linear form on (a subspace of) the linear space of bounded measurable functions 
ipn '. K”' —M with compact support defined by (we abbreviate “p.d.” for “pairwise distinct”) 

(2.4) Pn,p{,Tri) — Ep ^ ^ ¥^n(^l) • • • ; ^n)- 

Xl,...,Xn£C p.d. 

If the re-point correlation function exists as a distribution and can be identified with a measurable 
function, we will write f Pn,PTn instead of Pn,p{Tn)- Heuristically speaking, the one-point 
correlation function pi (also called the intensity of the random point process) gives the density 
of the process at each point, while the two-point correlation function p 2 {x, y) gives the probability 
of having a point both at x and y. In this paper we will work with stationary random point 
processes such that pi = 1 and we denote by Vsp{X) this set. 

2.4 Dimension extension 

We recall some elements from |PS15j to which we refer for more details. Outside of the Coulomb 
cases the Riesz kernel g is not the convolution kernel of a local operator, but rather of a fractional 
Laplacian. It can be transformed into a local but inhomogeneous operator of the form div(|ypV-) 
by adding one space variable y G M to the space In what follows, k will denote the dimension 
extension. We take /c = 0 in all the Coulomb cases, i.e. s = d — 2 and d > 3 or (IL2D and in all 
other cases we take k = 1. We use an auxiliary parameter 7 defined by 

(2.5) 'y:=s — d + 2 — k 

where the convention is to take s = 0 in the logarithmic cases. In particular we have 7 = 0 in 
the logarithmic cases (II.1|) (where k = 1) and (II.2p (where A: = 0). 
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Points in the space will be denoted by x, and points in the extended space by 

X = {x,y), X G y G The interaction kernel g is naturally extended to We 

will often identify x {0} and The measure S^d is the Radon measure on which 

corresponds to the Lebesgue measure on the hypersurface C 

Finally we let Cd,s be the constant depending on d,s, such that —div(|ypVg') = Cd,s^o in 
l^d+fc values of Cd^s nre given in |PS15j l. 

2.5 Electric fields and random electric fields 

Let p < be fixed. We think of the space as the space of electric fields 

generated by the charged particles together with a certain uniformly charged background and 
we endow this space with the weak topology. 

Local electric fields. If C is a finite point configuration and i? > 0 we let 

( 2 . 6 ) := Cd,s9 * (C - := = Cd,s^g * (C - IcAd). 

where g* denotes the convolution (computed in with the interaction kernel g. It implies 

that 

(2.7) - div(lyrEi-) = Cd,sVg * {C - IcAd). 

The scalar field physically corresponds to the electrostatic potential generated by the 
point charges of C together with a background of density g. The vector field can be thought 
of as the associated electrostatic field. It is easy to see that fails to be in because it 
blow ups like near each point of C, however is in . 

Electric fields. We now introduce a special class of vector fields that correspond to electric 
fields generated by a system made of an infinite point configuration C and a negatively charged 
background in all We let Elec be the class of “electric vector fields” i.e. the set of vector 
helds E belonging to M'^+^) that satisfy 

(2.8) - div(|ypE) = Cd,s {C - d^d) in 1^"^+^ 

where C £ X (M'^) is a point conhguration. We say that E is compatible with C if (12.81) holds. If 
E is in Elec we let 

Conf(E) := —3^div(|yPE) + 6^d 
C-d,s 

be the underlying point configuration C, in other words E is compatible with C if and only if 
Conf(E) = C. 

Truncation procedure. The renormalization procedure of |RS15j . |PS15j (inspired by the 
original work of |BBH94j 1 uses a truncation of the singularities which we now recall. We define 
the truncated Riesz (or Coulomb, or logarithmic) kernel as follows: for 1 > ry > 0 and X G 
let 

(2.9) UX) = [g{X)-g{g))^ 

with a slight abuse of notation: since s' is a radial function we write g{g) for the value of g at 
any point on a sphere of radius ry. 
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If is a local field as in (|2.6I) we let 

(2.10) := E^°^{X) - Y, V/,(X - p). 

p€C 

Similarly if E is an electric field as in (I2.8h we let 

(2.11) Er^iX) := E{X) - Y ^fviX - p)- 

p€C 

Random electric process. A probability measure on V concentrated on 
Elec is called a random electric process. We say that is stationary when it is invariant 

under the (push-forward by) translations E i—)• E{- — x) for any x G 

2.6 Specific relative entropy 

Let P be a stationary random point process on The relative specihc entropy ent[P|n] of P 
with respect to 11, the law of the Poisson point process of uniform intensity 1, is dehned by 

(2.12) ent[P|n] := , 

where P\Cii denotes the random point process induced in Cr, and Ent(-|-) denotes the usual 
relative entropy (or Kullbak-Leibler divergence) of two probability measures defined on the 
same probability space. We take the appropriate sign convention for the entropy so that it is 
non-negative: if p, v are two probability measures defined on the same space we let Ent (/u|z^) := 
Jlog if p is absolutely continuous with respect to u and -|-oo otherwise. We have in fact 
by super-additivity 

(2.13) ent[P|n] = sup Tv^Ent (P|c«|n|c^) . 

R>1 

The functional P i—)■ ent[P|n] is affine lower semi-continous on Vs^i(X) and its sub-level sets are 
compact. We refer to [RASOQl Chap. 6] for a proof of these statements. 

3 Definitions for the energy of a random point process 

In this section we recall the derivation of a renormalized energy for random point processes from 
(|1.4I) then we introduce our alternative definition of a (logarithmic, Coulomb or Riesz) energy 
for random point processes. 

3.1 The electric approach 

Renormalized energy of an electric field. We now recall the computation of the renor¬ 
malized energy W{E) following |PS15j (see also |Serl4] and the references therein). For any 
E G Elec we define 

(3.1) yVrjiE) := limsup (—i f l2/r|E^P - gip)) , 

il—>00 \Cd,s it / 

where is the truncated field as in (12.1111 . and we let 

(3.2) >V(P) := limWJP). 

For convenience we have chosen a definition of W which differs from a multiplicative constant 
Cd^s from the one in |PS15j . 
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The electric definition. We then define 

(3.3) W^(C) = inf>V^(S), W(C) = inf >V(.E), 

where both infimum are among electric fields E compatible with C. Similarly if P is a random 
point process we let 

(3.4) W®^®'=(P) = Ep [W^], W®^®'=(P) = Ep [W]. 

The following lemma was proven in [LS15j : 

Lemma 3.1. Let P be a stationary random point process such that W®^®'^(P) is finite. Then 
there exists a stationary random electric process such that the push-forward of P®^®® by 

Conf is equal to P and which satisfies 

(3.5) Epelec[W] = W®^'^'=(P). 

We also have the following lower semi-continuity result for the electric energy. Let us em¬ 
phasize that lower-semi continuity only holds at the level of stationary random point processes, 
and not for point configurations or arbitrary random point processes. 

Lemma 3.2. The maps P i—)• W®*®'^(P) and P i—)• W®^®'^(P) are lower semi-continuous on 

Proof. See e.g. |PS15[ Lemma 4.1.] □ 

This provides a definition for the energy of a random infinite point configuration with a 
uniform negative background. However the computation of W*^^®‘^(P) or of W(C) (or even the 
search for an upper bound on these quantities) appears involved in general because it amounts 
to finding compatible electric fields for infinite point configurations. Let us mention that in the 
case of a periodic configuration (hence also for the stationary random point process associated 
to it) exact formulas are known in the cases (11.21) (see Esm), 0 (see |SS14] i. for the higher¬ 
dimensional Coulomb case m with s = d — 2 (see |RS15j i and for Riesz gases (see |PS15] ). 


3.2 The intrinsic approach 

In this section we define an energy functional W™* on the space of random point processes, using 
only the nature of the pairwise interaction. 

If A,B are two (measurable) subsets of we define Int[74,R] as the interaction energy 
between A and B 


(3.6) 


lnt[A,B]{C) 


{AxB)\A 


g{x 


y){dC{x) - dx) (g) {dC{y) - dy). 


In view of (m, for any R > 0 we define TTf^ : X{Cr) —)• M as the interaction of Cr with 
itself (the diagonal being excluded) i.e. TTf^ := Int[C'ij, C/j] or in other terms 


(3.7) Wr{C):= // g{x — y){dC — dx) ® {dC — dy). 

JJcj^\A 

Given a random point process P, a natural way of defining the energy (per unit volume) of P is 
the following: 

Definition 3.3. Let P be a random point process of intensity 1. We define its intrinsic energy 
by 

(3.8) W”*(P) := Immf ^Ep [n^\C)] . 
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Expression with correlation functions. The energy defined by (13.81) can be re-written with 
the help of one- and two-point correlation functions of P. In order for the expression to make 
sense, we restrict ourselves to random point processes whose two-point correlation function exists 
as a Radon measure in x We will then abuse notation and consider p 2 as a function 
instead of a measure, writing p 2 {x,y) instead of dp 2 . 

Lemma 3.4. For any random point process P of intensity 1 such that p 2 exists as a Radon 
measure, the following identity holds 

(3.9) W“*(P) = liminf ^ /7 g{x - y){p 2 {x,y) - l)dxdy. 

R^oo JJc%\A 

Proof. First we may re-write TPf^{C) as : 


'1/int 


(C) = // g(x — y)(dC — dx) (8) {dC — dy) 

JJcl\A 


g{x — y){dC (8) dC) -|- 


g{x — y){dx ( 8 > dy) — 2 


c%\A JJcl\A JJcl\A 

and by definition of the correlation functions (see Section [2.3p we get 


g{x — y)dC (8> dy 


Ep 


['H'r]= g{x-y)p 2 {x,y)dxdy+ g{x - y)dxdy - 2 g{x - y)pi{x)dxdy. 

JJcl\A JJcl\A JJci\A 


<Cl\A J.Jcl\A JJcl\A 

By assumption P has intensity 1 i.e. pi = 1 and we are left with 


[Wk“] = fs 


R'^ ././C|\A 


9{x - y){p 2 {x,y) - l)dxdy, 


which yields (13.9p . 


□ 


The stationary case. If the random point process P is stationary we may derive a somewhat 
simpler expression for W“*(P). In what follows P 2 {v) stands for /92(0, v). The change of variables 
(tt, v) = {x + y,x — y) gives 


{l-R/ 2 ,R/ 2 ]‘i)‘^\A 


g{x-y){p 2 {x-y)-l)dxdy =-^ [ [ g{v){p 2 iv) - l)dudv 

2 JvG\—R.R]‘^\tO\ J u£Sr(v) 


j;G[— JueSii(v} 

\SR{v)\g{v){p 2 {v) - l)dv, 


I 


vel-R,R]d-\{0} 


where Sr{v) denotes the set 

Sr{v) := {x + y : x,y £ Cr, X - y = v}. 

The Lebesgue measure IS’pPu)! of Sr{v) is easily computed for v = (ui,... ,Vd) 


(3.10) 


|5p(ui,...,urf)| =2''x[](R-|ui|). 


i=l 


Indeed we have Sr{v) = {2x — v : x,y £ Cr, x — y = v} which implies that 

\Sr{v)\ = 2 X |{x G Cr,x- V £ Cp}! , 


13 


moreover Sji{v) tensorizes i.e. if ?; = (?;i,..., Vd) G [—-R, we get 

d 

(3.11) e [-i?/2, ii/2], x, - e [-R/2,R/2]]\ 

i=l 

which leads to (I3.10p . Finally we obtain the following expression for any P in Vs^i{X) (such 
that p 2 is a Radon measure) 


(3.12) 


wmt(p) ^ 1 f g{v){p 2 {v) - 1) f[{R - \Vi\)dv. 

R^oo R J[-R^R]d\{0} 


4 Preliminary results on the energy 

4.1 Local field and local interaction 

Let R ^ 0 and let d be a point configuration in ^^(Cr^. The local electric potential (resp. field) 
(resp. are defined in (j2.6p . For p S (0,1) we let also be as in (12.10^ . Let us recall 
that Cr denotes the hypercube [—R/2,R/2]'^^^. 

Lemma 4.1. The following inequalities hold 

1. In the cases don 

(4.1) — [ \yr\E^^^\^-MRg{p) < nf[C] +Mro,{1). 

^d,s 

2. In the logarithmic cases 


(4.2) —[ \y\^\E^^^\‘^-MRg{v)<nf[C] + &°^VllogR + MRO,il) 

C-d,s JCu 

+ 0(MrR ^) + Oi?(l). 


Let P be a stationary random point process of intensity 1 such that Ep[P^] = o{R?'^). The 
following inequalities hold 
1. In the cases (IL31) 


(4.3) 


Ep 


1 

C-d,s 



■f^Rdiv) 


<Ep[nf] + R\{i). 


2. In the logarithmic cases 


(4.4) Ep 


1 


— I \y\^\E^f^r-^fR9iv) 

('d,s JCn 


< Ep[n^^] + &^^Ep[Vl] logR + R'^Orjil) 

+ op(l). 


The terms o^(l),op(l) depends only on d,s. 

Proof. The starting point is the following identity which holds for any S > R 

(4.5) W^^[C] = lim f— [ \y\^\E\°^\^-MRg{p)) + — [ • n, 

^^0 \Cd,s Jcs ) Cd,s JdBs 
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where n denotes the unit normal vector. It results from two different operations. Let us first 
recall that 

•= ff g{x — y){dC — dx) ® {dC — dy) 

JJcl\A 

Since satishes (12.7h we may formally write 

'Hr"[C] ~-^ f d>'°'^(t)div(|yrV$'°")(t). 

Of course div(|y|'^V<h^°^)(t) is in fact supported on Cr. For any 5 > i? an integration by parts 
yields 

nf[C] ! | 2 /r|L;'°"P + — • n. 

C-d,s J Cg C-d^s J dCs 

Since is not in the previous computation does not make sense, however it can be made 
rigorous in a renormalized fashion by truncating the interaction close to the charges at scale y > 0 
and substracting a diverging term as r/ —^ 0 (here g{r])) for each charge. We refer to [PS 15] for 
more details. 

We now turn to the boundary term in (14.51) . A mean value argument applied to g and Vg 
shows that for S > 2R we have 


_ VRgiX) < C{Nr + R‘^)RS 


— S—1 


\E^--[X) - VRXg{X)\ < C{Mr + R<^)RS-^-\ 

uniformly for X S dCg- Indeed it is easy to see that the first derivative of t i—>• g{X — t) is 
bounded by CS~^~^ and that its second derivative is bounded by CS~^~‘^ for any t in Cr, 
uniformly for X G dCs, and with a constant C depending only on d, s (because S > 2R). 
Moreover we have |V( 7 (S')| < the perimeter of dCg is 0(5”^+^“^) and \y\^ < CS'^ on 

dCs- Since 7 = s — d + 2 — A:we have 5'7+rf+*:-i = 

Combining the estimates above we get for any S > 2R 



|y|7$loc^loc 


n 


< ( V 


'R' 


| g(^)l 

5 .+I 


+ (A/1 


1 + -^ 


2d\ 


R , ^|g(^)l 

^ 2 s +2 5«+2 


—'ll 

S2.+3 J ) 


with a constant C depending only on d, s. 

In the cases (|1.3p the right-hand side is 0(5“®) as 5 ^ 00 and s > 0. We thus obtain 


(4.6) 


[ • n 

Jacs 


— 05_>.oo(l)- 


In the logarithmic cases we have s = 0 and g{S) = log 5, hence if we set S = R'^ we get 

(4.7) [ • n < logi? + O {Af^R-^) + oh(1), 

JdC^4 

for a certain universal constant 0*°®. 

Combining (|4.5p and (14.61) and letting 5 —)• 00 we obtain, in the cases ()1.3p 


lim(— [ |!,n<'P-VM(i,)) =Hg‘[C|- 

7^0 \Cd,s jRd+k J 
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Combining (14.51) and (14.71) we obtain, in the logarithmic cases 

hm f— [ \yP\E^°^\^ -Mngiv)) < [C] + &°^Vllog R + O + or{1), 

^ 7^0 \ Cd,s JCr ) 

where we have used the trival bound Ini'll^ \yy\E^°^\^. 

From |PS15[ Lemma 2.3] we know that the limit as r/ —)• 0 is almost monotonous, more 
precisely for any rj G (0,1) and S > 2R we have 

(4.8) [ \y\^\E^°^\‘^-AfRg{y) < lim (— f \y\^\E';;‘^\^ - AfRg{g)\ +AA«o,(l). 

J(Jg ^ ^ JCq J 

We thus obtain (HU) and (1121). Inequalities (1131) and (1131) follow easily by taking the 
expectation under P. □ 


4.2 Bound below on the interaction 

The self-interaction of an electric system of point charges with a negative background in a given 
set K is bounded below in terms of the number of points in K (and the volume of K in the 
logarithmic cases). 

Lemma 4.2. Let K be a compact subset such that K C Cr and let C be a point configuration 
in K. We denote by M{C,K) the number of points in K, and we denote by Int[iL, iL](C) the 
interaction energy of K with itself as in 
In the cases dEH) we have 

(4.9) Int[iL, K] (C) > -C {J\f{C, K)) . 


In the logarithmic cases we have 

lnt[K,K]{C) > -C[M{C,K) + {Af{C,K) - \K\)^ log R + Af^{C, K)R-^) . 


In both inequalities C is a positive constant depending on d,s. 


Proof. Letting E^^^^ be the local electric field generated by C in iL and arguing as in the proof 
of Lemma O we obtain 


lnt[K,K]{C) = ]im(— [ \y\fiE];>^\^ - AfRg{g)) +— [ ■ n. 

\Cd,s Jcs J Cd,s JdBs 

As seen in (14.81) the limit as r/ —)• 0 is almost monotonous. In the cases (II.3p . chosing g = ^ 
and letting 5 —)• oo we obtain 

(4.10) lnt[K,K]>— f \yP\E^n^-M{C,K){g{l/2) + C), 


which yields (14.9p . 

In the logarithmic cases, choosing g = ^ and taking S = R‘^ we obtain 


/ |!/r|M“p-V(C,/r)(9(l) + C)- 


/c,4 » 2 

JdC^4 2 


Controlling the boundary term as in the proof of Lemma l4.II yields pd.lOp . 


□ 
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4.3 Discrepancy estimates 

The following lemma is an adaptation of the discrepancy estimates of [PS15j to show how the 
hniteness of W™*(P) (instead of implies that P has a number variance of order 


Lemma 4.3. Let P be a stationary random point process of intensity 1. 
In the cases do]) ifW™^{P) is finite then we have 

(4.11) Ep[Vl] = 


In the logarithmic cases ifTf°^{P) is finite then we also have (14.111) . 

Proof. In the logarithmic cases, the result is straightforward since T>^°^{P) < +oo implies that 
Ep[Vl]=o{R<^). 

We now turn to the cases dO]). Applying |PS151 Lemma 2.2] with jy = ^ we get 


P 


:V%mm ( 1, 


pR 

R<i 


Taking the expectation under P yields 


Ep 


1 


T>P min ( 1 




< C 


I 

Jm 


Kd+fc 2 


iloc 1 2 


+ CMr. 


< CEp 


jRd+i; 2 


loc |2 


+ CR‘^, 


where denotes the local held generated by C in Cr and C is a constant depending on d, s. 
We know from Lemma l4. II that 



J\fRg{\)<R}R^[C] + CJ\fR, 


and combining the previous two estimates we obtain that 


E, 


1 -n2 ■ f. Pr^ 


R<i 


< CEp[nf] + CR‘^ = 0{R^ 


Using Jensen’s inequality and the fact that s < d we get 

Ep[p|] = 0(iit'^+i") = o{R^'^), 


which proves ()4.1ip . 


□ 


4.4 The screening lemma 

For convenience we recall the “screening lemma” following [SS12] . |SS14j . [RSI 5] . |PS15j and 
[Lsn]. We present it here in a simplihed form which will be enough for our purposes, we refer 
to [LSI 5] for the most general statement available and to [PS 15] for a proof. 

The result consists in the following: given a point conhguration C in Cr and a compatible 
electric held E, we wish to construct another compatible held E^^^ such that • n = 0 on 
the boundary of Cr x M^. Indeed such screened helds may be pasted together in adjacent 
hypercubes because their normal component vanish (and in particular are equal), which will 
allow us to construct global helds dehned in the whole space. 

Of course this construction is not possible in general (e.g. it is easy to see that it imposes a 
condition on the number of points in Cr, which must match exactly the volume of Cr). However, 
under some conditions on E to be “screenable”, by extending Cr a bit and modihying C only in 
a thin layer of width eR we may hnd a new point conhguration C®™ and a compatible screened 
held E^^'^, such that moreover the energy of E^^'^ is bounded in terms of the energy of E. 
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Lemma 4.4. There exists Rq > t) depending on d, s and r]Q > 0 depending only on d such that 
the following holds. 

Let 0 < e < ^ and 0 < rj < rjQ be fixed. Let Cr he a hypercube of sidelength R for some 
R > 0, and let K be the hypercube of sidelength [i?] (where [i?] denotes the smallest integer 
larger than R). 

Assume that E is a vector field defined in Cr x Mf such that 


-div(|yP.E) = Cd,s {C - 
Let M > 1 such that E satisfies: 


in Cr. 


IcrxI 


\y\fiEr,\^<M. 


In the case k = 1 (the non-Coulomb cases) we define e^^R as 


(4.12) 


ee,R ■= ^4 




\y\^\E\^- 


ICRX(S.\{-e^R,e^R)) 

Under the assumption that the following inequalities are satisfied 


(4.13) 


R > max 


Rq RqM \ 


R > 


RoM^/^ 

gd+3/2 


2 j_ 

,RQe^-~<e 


l/(l-7) 

e,R 


if k = 0 


ifk = l, 


there exists a point configuration in K and a vector field E^^^^ G M'^'*'^) such that 

1. The configuration has exactly \K\ points in K. 

2. The configurations C and coincide on Intg := {x G C/j, dist(x, SC/j)} > 2ei2}. 

3. We have 


I -div (|yr£;^“) = Cd,s (^^“ - 4^) tnKxR^ 
l^scr.-^O ondKxR^, 


4 . Letting be associated to as in (12.111) it holds 


(4.15) [ \y\^\E)(r<([ |yn.L,|2)(l + C6) 

JKxR>= V./CflXRfc / 

+ Cg{r])MsR‘^ + Cee^RsR'^ + o{R!^-^). 

for some constant C depending only on s, d. 


4.5 Minimality of the local energy 


The following was proven in [LSlSi Lemma 3.13.]. It expresses the fact that the local electric field 
has a lower energy than any “screened” field compatible with the same point configuration. 

Lemma 4.5. Let R > 1 be an integer and C be a point configuration in Cr. Let be the 
local electric field generated by C in Cr as in (12.61) . Let E G ,R^~^^) be a vector field 

satisfying 


( -div(|y|'^£;) = Cd,s {C - (^Rd) in Cr x R^ 
\ E-11=0 on OCr X 


Then, for any 0 < rj < 1 we have 


(4.17) 



< 



2 
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5 Connection of the electric and intrinsic approach 

This section is devoted to the proof of Theorem [TJ It goes in two steps. 

1. First we establish an upper bound < W™* (or W™* + T>'°s in the logarithmic cases). 

The proof of this “electric-intrinsic” inequality is the purpose of Section 15.11 

2. Then, in the non-Coulomb cases, for any P S Ps,i{^) we construct a sequence {Pn}n 
of stationary random point processes which converges to P and such that W®^®®(P) is 
bounded below by lim^v-^-oo W™*(PAr) (we also ensure that Pn is hyperuniform, in partic¬ 
ular P^°®(P/v) is always zero). Moreover we have limAr_>.ooent[P 7 v|n] = ent[P|n]. 

This operation is similar to the construction of a “recovery sequence” in T-convergence, 
and is proven in Section 15.21 
These two steps immediatly imply Theorem [TJ 


Re-writing of the additional term. Using the definition of the two-point correlation func¬ 
tion (which exists by assumption as a Radon measure) we get 


Cl 


{P 2 {x, y) - l)dxdy = Ep [AAp(AAp - 1)] - = Ep [AA|] - - Ep [AAp]. 

Since P has intensity 1 we have Ep [A/p] = R'^ hence 

(5.1) [[ ^iP 2 {x,y) - l)dxdy = Ep [P|] - R‘^. 

Equation (|5.1I) allows us to write the term (defined in (ll.Sp l in terms of p 2 equivalently as 

(5.2) P'°s(P) = (^'“Shmsup I ^ f[{p 2 {x,y)-l)dxdy + l\logR 

R ^oo \ JJ j 

= lim sup ^Ep [p| log R] . 
R-^cx> ^ 


5.1 The electric-intrinsic inequality 

Until the end of Section 15.11 P denotes a stationary random point process of intensity 1 on 
such that W™** is finite. In the logarithmic cases we assume that P^°^{P) is finite. 

The following proposition is the hrst part of the proof of Theorem [TJ 

Proposition 5.1. Under the above assumptions, we have in the cases m 

Wclec(p) < 

and in the logarithmic cases W®*®®(P) < W™*(P) -T 'D^°^{P). 

Screenability of the local electric fields. Let ? 7 , e > 0 be fixed. Let {Rn}n be an increasing 
sequence of real numbers such that lim„_).oo Rn = +oo and 

(5.3) lim ^Ep[P”*]=W“*(P). 

n^oo Kf, 

We start by an auxiliary lemma. 
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Lemma 5.2. Let denotes the local electric field generated by a point configuration in Cr^. 

The following inequality holds 


(5.4) 


P 


L 




|y|7|^loc,n|2 < ^ ^ ^ _ 0 {M-^), 


moreover, in the case k = 1 , 
1 


(5.5) P 


L 






27 \ 7-1 


R 


d—s 


= 1 - o„(l). 


Proof. The first point (|5.4p follows directly from Markov’s inequality and the assumptions on 
P. Indeed we have, in the cases (11.31) (using (I4..sp l 


Ep 


— f 

Cd,s JCfl„xR'= 


<Ep|w£] + 0(j;“), 


and in the logarithmic cases (using (14.41) 1 


Ep 


- f i2/riEf 

Cd,s Jc'h„xK'= 


< Ep[ll^^:] + Ci°gEp[P|j logi?,, + OiRi), 


where the terms 0 {R'^) depend on g. 

To prove the second point (15.5p let us fix X = {x,y) with x G Cr^ and y > e^Rn- We 
may estimate E*°'’’"'(X) as follows: let i?o > 0 and let us divide Rn into 0 {R'^/Rq) hypercubes 
{Ci}i^j of sidelength G {hRo, |.Ro)- For any i G / we have by a mean value argument 


[ Vg{X - t){dC{t) - dt) 

JCi 


<c\v. 


l-s-1 


+ CRo{Mi + Rt)\y\ 


d\\^.\—s—2 


where Pi (resp. A/)) denotes the discrepancy (resp. the number of points) in Ci. Summing over 
i £ I we get 

|^loc,n(^)| < ^ CRo{Mr„ + R^I.M- 

iei 

We may thus write, using the stationarity of P 


\-s-2 


M: 


2d 


2d\\„,\-2s-A 


Ep[|E'-’-(X)|2] < C^Ep[P|j|y|-2-2 + Ci?o(E[AA2j + Rl^)\y\ 

Rq 

Using Lemma [4.111 or the fact that ^^“^(P) is finite we have E[A/p^] = R^ + 0{R^). Finally, 
integrating over Cr x M\{—e^Rn,£^Rn)) we obtain 


Ep 




■n ./Cij„x(K\(-e 2 p„,£ 2 p^)) 


|y|7|^loc,n|2 


< ^ CEp[P4 


T32d 


' Rf {e‘^RnY+^ 


+ C RqR, 


2d 


(e2p„)'^+^+2y ■ 


Using the fact that ■j^'Ep[T>‘j^J = o(l) as Pq oo (see Lemma IT^ . we obtain 


E p 


'Rq 


L 


^‘^Rfi x(R\(-£ 2 p„,e 2 p^)) 


|y|7|Eloc,n|2 


= o{R. 




The bound (15.51) follows by Markov’s inequality. 


□ 
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We now turn to the proof of Proposition 15.11 


Proof. Screening the local electric fields. When r],£ > 0 and M > 0 are fixed we denote 
by the set of point configurations in such that satisfies (15.4p and, in the case 

k = I, (15.51) . For any C in the conclusions of Lemma 14.41 apply to . We may thus 

find a point configuration in := and a compatible field E^^'^ such that 

1. The point configurations C and coincide on a large subset of Cr^., namely {x G 
Cr„, dist{x,dCR„ > 2£Rn}. 

2. The vector field is screened i.e. E^^^ • n = 0 on dKn x M^. 

3. The energy of E^^'^ is bounded in terms of that of as in (|4.15p 

(5.6) [ <([ (i+ 

+ Cg{v)MeRi + Cee,R£Ri + oiR^), 


where e^^R is defined in (14.121) . 

Constructing a global electric field. Any such point configuration (resp. vector field) 
may be extended periodically in the whole space The main point is that since E^^^^ is 
screened we can paste together several copies of E^^^ periodically without creating divergence 
at the boundary of two tiles. Let CP®’’ (resp. F^p®’’) be the resulting periodic point configuration 
(resp. vector field) 

We let Pn be the conditional expectation of P knowing and we let be the push- 

forward of Pn by the map C i—^ PP®’’ defined above. Finally we let Pf^^ be the average of Pn^^ 
over translations in Kn- Taking the expectation of (j5.6p under P we see that (with W®'®® as 
defined in (pS.dp i 


w®^®®(pr) < 


<C -tEp 

- Ri " 


-f 

Cd,S J KnXi 

-/ 

i,s J C 




-9{r]) 


Cd,s JCr^xI 




(1 -h Ce) - g{r]) + Cg{g)M£ 

+ CeBp^les^Rj -I- o{Rn^). 


From Lemma 15.21 we see that as M ^ oo, n —)• oo the random point process Pn converges 
to P. In particular we may bound the expectations under Pn in the right-hand side by the 
expectation under P at a small cost 



(1 + 0 ( 1 )), 


< Ep[e£,K](l + o(l)) 

where both terms are o(l) as M —>■ oo,n —>■ oo (keeping e,g fixed). By definition of and a 
mean value argument we see that up to changing e into 2e we may assume that 


Ep[e, 


£,Rn 


= Ep 






|yriE 


iloc,n|2 


< 


e6p, 


-M. 
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We thus get 


wf'=(pr) < 


-f 

C-d,s Jc 






(l + o(l))(l + Ce)- 5 r(r/) 

+ Cg{r])M£ + C£~^R~^. 


The convergence of to P implies the convergence of to P as M, n ^ oo and e ^ 0. 
This might be seen as follows: the topology on V{X) is such that for any fixed 5 > 0, if two 
random point processes coincide (or are close to each other) in Crq for Rq large enough, then 
they are 5-close. In particular Pn and P are very close to each other in Cr^ because on the one 
hand the vast majority of point configurations under P are screenable and on the other hand 
the screening procedure does not modify the points in a large interior part of Cr^ . Heuristically 
speaking, if Rn is larger than some Rq then and P should be 5-close. Now when averaging 
the random point process Pn over translations in Cr^ there is a nonzero proportion of z ^ Cn 
which are such that the translation by z of the thin layer of Cp„ in which the points have been 
modified ends up intersecting Crq, thus P^'' look less like P in Cq. However when Rn S> Rq 
this proportion is of order e, hence P^'' still converges to P when taking M, n large and e small. 

Conclusion. Taking M large, letting n —)• oo and using (14..SI) or (14.41) we obtain, by lower 
semi-continuity of over Vs^i{X), (see Lemma 15?^ that 

W®i®®(P) < W”*(P)(1 + Ce) + 0(e) + o^(l), 

plus an additional (D^°®(P) term in the logarithmic cases. Sending e ^ 0 and ?] —)• 0 we finally 
obtain that W®^®®(P) < W™*(P) or, in the logarithmic cases W®^®®(P) < W“*(P) -|- P^°®(P), 
which concludes the proof of Proposition 15.11 □ 


5.2 Construction of a recovery sequence 

In this section P denotes a stationary random point process of intensity 1 such that W®^®®(P) is 
finite. Moreover we assume that we are in one of the non-Coulomb cases, i.e. d = 1 or d >2 
and s > d — 2. 

The following result forms the second step in the proof of Theorem [TJ 

Proposition 5.3. There exists a sequence {Pn}n of hyperuniform random point processes in 
Ps^i{X) such that 

lim P]\f = P, lim ent[P 7 v|n] = ent[P|n], 

N^OO N—^OO 

and satisfying 

(5.7) lim W“*(P 7 v) = W®^®®(P). 

TV —>-00 

Let us observe that since the random point processes are hyperuniform, in the logarithmic cases 
they satisfy V^°^{P]\f) = 0 for any N. 

The proof of Proposition 15.SI goes in two steps. 

1. First we construct an auxiliary sequence of random point processes which converges to P 
and such that almost every point configuration is finite and “screened” i.e. there exists an 
associated screened electric field. This is done in Lemma 15.41 

2. Next, we extend this random point process in the whole space and make it stationary, 
before bounding its interaction energy in terms of W®^®®(P). 
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5.2.1 An auxiliary sequence 

Lemma 5.4. There exists a sequenee {P^^}n of random point processes in V{X{Cn)) sueh 
that 

0. The sequence {P^^}n eonverges to P as N ^ oo. More preeisely, there exists a sequence 
{L]y}]y sueh that Ljy = N{1 — o(l)) and the respective restrictions of and P to Clj^ are 

arbitrarily close as N oo. 

1. For P^'^-a.e. point eonfiguration there exists a screened eleetrie field satisfying 
(5.8) - in Cn x 


(5.9) • re = 0 on SCn x M^. 

In partieular the point configurations have P^^-a.s. N'^ points in Cn- We also have 

(5.10) min dist(p, SC^v) > rjQ, 

pGC(l) 


for some % > 0 depending only on d,s. 
2. The following estimate holds 


(5.11) 

3. 

(5.12) 


lim sup lim E 


(1) , I 

N^oo lCd,s Jcj^x 


1 1 


[ \yr\Pv^\‘^ - aiv) 

JCm xM^ 




The relative entropies of P^^ and P\Cn with respect to are elose 


Ent[Pi')|n|c,] = Ent[/)c7.|n|c;,] +o(A''). 


Proof. This follows from the analysis of |LS15j . and we sketch here the main steps. 

Let be a stationary electric process associated to P as in Lemma 10.11 For fixed 

R, M, e, r/ > 0 we say that an electric field E is in (or is screenable) if its energy is 

controlled as follows 


1 

'W 



< M and, if A; = 1, 




L 


CRxRfi(-e^R,e'^R) 




< 1. 


Under the assumption that W®^®'^(P) is finite, then the probability tends to 1 as 

M, R ^ oo for any e, ry > 0 hxed. This is proven in |LS15l Lemma 5.10] and is similar in spirit 
to Lemma 15.21 

If R is an integer much larger than M and E is in , the screening procedure as in 
Lemma [13] (see also |PS151 Proposition 6.1] and |LS151 Proposition 5.2]) applies. In particular 
we may change the underlying point configuration in a thin layer of size < eR close to the 
boundary of Cr and obtain a new sereened point configuration C®®'’ in X{Cr) as well as a 
compatible screened electric field E®®'’ which satisfies (15.9p . It also ensures that (I5.10p holds. 
The screening procedure is described in [LS151 Section 5.1]. The energy of E®®’’ is bounded in 
terms of the energy of E as in (I4.15p . 

The next step is to regularize the point configurations C®®'" by separating the pair of points 
which are close from each other. This regularization procedure is described in |LS151 Section 5.2], 
and another electric field E™°®^ can be associated to the regularized point configurations, with 
a good energy bound. The main benefit of this procedure is to control the difference between 
W /cflxRfe - aia) and its limit as ?7 —)• 0. In general the limit may be much larger 
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because of the contribution of pairs of points which are very close, at distance <C rj, and which 
are not “seen” when truncating at scale r] > 0. 

We let e, r] tend to 0 and M tend to infinity (depending on N) and we pick R = N large 
enough. We let be the associated random point process in Cn- Most of the point configura¬ 
tions (or electric fields) are “screenable”, the screening procedure only modifies the configuration 
in a thin boundary layer of Cw, and the regularization moves only a small fraction of the points 
by a small distance. This ensures that converges to P as —)• oo. 

The estimates on the energy of the screened-then-regularized electric fields are such that 
(|5.1ip holds (see |LS151 Section 6.3.4.]) with := 

Concerning the entropy, letting the new/modified points of the configurations move in small 
balls allow us to recover a small, nonzero volume in phase space without affecting the energy, 
since only a small fraction of the points have been deleted/created/modified, it gives (I5.12p 
(see [LS151 Section 6.3.5] for a precise analysis of the volume loss). □ 


5.2.2 Proof of Proposition ISTHl 

Proof. Step 1. Construction of the random point process. Let be as in Lemma 15.41 and 

let us extend P^^ in the whole space as follows. Let be a tiling of by a family of 

hypercubes of sidelength N and let Xi be the center of Xi (we may impose that one of the xfs 
is 0). Let be the laws of independent point processes distributed as P^^^. 

To any family ig/ of point configurations in Cn we may associate the point configuration 

C:= 

iei 

which amounts to “paste” the point configuration in the hypercube Ci. 

For any i ^ I, let be an electric field which is compatible with as in (15.8]) and 

screened as in (15.91) . By the latter condition, the normal component of each vanishes on the 

boundary of Cjv x thus we may paste together such fields along the boundaries and their 
energy is additive. In particular the electric field P defined by 

(5.13) P(x) := ^pW(a;-Xi) 

iei 


is compatible with C and moreover we have 




for any measurable subset A C and any rj > 0. Let us also observe that, by construction, 
the normal component of P vanishes on the boundary of CmN x for any m > 1 (in fact it is 
easy to see that it vanishes on any path included in x M^). 

We let P^^ be the random point process obtained by pasting P^] in Ci for i S /, or in other 
terms the push-forward of the product measure of the by the map 

i&I 

(2) (2) (2) 

For any z G Cwj we let Pff'^ = Oz ■ Pff' be the push-forward of Pf^' by the translation by a 

( 3 ) ( 2 ) 

vector z. We also define P^ as the uniform average of Pfi'^ for 2 : G Cw- It is not hard to check 
( 3 ) 

that P^' is a hyperuniform stationary random point process which converges to P as N ^ 00 . 
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Step 2. Estimates on the energy. We fix z G Cat and m > 1, and we use the subscript z to 

denote a translation by z, e.g. CmN,z denotes the hypercube CmN translated by z. 

( 2 ) 

By construction, to ^-almost every point configuration in C^rn+i)N,z we may associate an 
electric field E whose normal component vanishes on the boundary of C(^m+i)N,z ^ 

Let us recall that Int[^, B]{C) denotes the interaction energy between the sets A and B (cf. 
(I3.6p l. By minimality of the local energy (cf. Lemma l4.5p we have for P^(,-almost every C 


lnt[C(^rn+i)N,z, C'(m+i)iV,2](C) < lirn ( — / \Er,\‘^ - N’^gi'n) 

Using (15.lip we see that the right-hand side is bounded in terms of as —>■ oo, more 

precisely we obtain 



1 


(5.14) limsuplimsup 

N—>-oo m—>-oo {mN) 




and both limits (as m —)• oo and N —>■ oo) are uniform for z £ C^. 

This gives an asymptotic upper bound on the expectation of Int[C'(m+i)Ar,z; C'(m-i-i)7V,2] un- 

der however the relevant quantity to control in order to get p 5 . 7 p is rather the expec¬ 

tation of Int[Cmiv, UmAr]. We thus need to bound the difference Int[C'(m+i)Ar,2) C'(m-i-i)Af,2] “ 
Int [CmN ) CmN] ■ 

Let us write Int[C'(m+i)Ar,2, C'(m+i)iv,2] as 


(5.15) ^'^C\C(m+l)N,z-:C(m+l)N,z\ — fat [UmAf j C*mAf] T 2Int [C*mAf; C'(m-|-l)Af, 2 \dmAf] 

-|- Int[C*(,.„^]^)Af,2\C*mAf) C(m+\)N,::\CmN]- 

We may bound the last term in the right-hand side of p5.15p using Lemma 14.21 with (for the 
notations of the lemma) K = C(^m+i)N,z\CmN and R = {m + 1)N. In C(^m+i)N,z\CmN there 
are points, and the discrepancy between the number of points and the volume is also 

such that {N{C,K) — \K\f‘ = 0(m'^“^). Applying (14.9p in the cases pi.3p and p4.10p in the 
logarithmic case we obtain 

(5.16) Int[C'(m+i)Ar, 2 \C'mAr, C'(m+i)iV, 2 \C'mAr] > o(m'^) as m oo. 

We are left to estimate the interaction term Int[C'mAr,U(m+i)Af, 2 \C'mAf]- We may split it as 


(5.17) RA\CmN ■iC(^m+l)N,z\CmN] fat [C7(^_]^)Af,2! C*(m+l)A'',2\C*mAf] 

-|- Int[C*.niAf\C*(m—1)W,2) f^(m-|-l)Af,2\f^mAf]) 


and we will prove 

(5.18) E (2) [lnt[C'mAr,C'(m+l)Ar,2\C'mAr]] = o(m'^). 

^N,z 


Studying Int[C'(^_i)Ar, 2 , C'(m+i)Ar, 2 \C'mAr]- First we claim that 
(5.19) E ( 2 ) [lat[U^^_]^)A/', 2 ! f^(m+i)Af, 2 \f^mAf]] oijn ). 

N,z 

To prove p5.19K let us write 

Int [C'(^-i)Ar,2,C'(m+i)iV,2\C’m7v] =- — [ $“div(|ypE^), 

Cd,s 
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where is the local electric potential generated by the system of charges in C'(m+i) 7 v, 2 \C'mAf 

i.e. 

:= f g{x — t){dC{t) — dt) 

(-'(m+l)N,z\CmN 

and E’^ is the screened electric field associated to the system of charges in C(^m-i)N,z- Using 
(|5.10p we may also write 


1 




Cd,S J C^jri-l)N,z 


ch^divdyruf 


970/) 


because the minimal distance between a point charge in C(^rn-i)N,z one in C(^ra+i)N,z\CmN 
is > An integration by parts and Cauchy-Schwarz’s inequality yield (observing also that 
C(m—l)N,z U m +2 )n) 


(5.20) Int[C(^_i)7v,z,C 

(m+l)JV,2 \CmAr] 

> -c 


1/2 


' C'(m+ 2 )ivXl 


mEi \^ 


fRd+fc 




1/2 


( 2 ) 

for some constant C depending on d, s. Using the definition of Plj ’ we have 


Ep(2) 

-* \r - 




\e: 


~b |2 
VO I 


< m'^Ep(2) 


where jg as 


m 


Rd+fc 

Using (I5.1ip we obtain that 
1 


/ 

JCm xR^ 


limsup E ( 2 ) 


\e: 


b |2 


A'"— 7 oo,m—>oo {jnN^^ J^d+k ^ 

with C depending only on d,s. In particular this yields 


< (W®^®^(P) + C) 


Ep(2) 


/ 

Jr‘*+'= 


\e: 


b |2 
VO I 


= 0(m'^). 


Using Jensen’s inequality we may thus bound the second term in the right-hand side of (j5.20l) 
as 


(5.21) 


E 


7| pnft 12 

VO I 


\ym: 


1/2- 




It remains to control the the first term in the right-hand side of (|5.20p . We claim that 


(5.22) 


Ep(2) 

N,z 


JC'(m+2)JV Xl 


\yV\E, 


VO I 


= o{m’^). 


Since is the local field generated by the configuration in C'(m-i-i)Af,z\UmAf we may, by almost 
monotonicty as in Lemma l4.II compare + 2 )jvxR'= with the interaction energy 


Int [C'(ni-|-l)Af,z\U'mAf) ) 
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up to a boundary term in the logarithmic case (which is bounded as usual) and a term of order 
0{m'^~^)g{r]Q) = o{m'^). We now prove that 

(5.23) E pi2) [^I^t \^^7?2A^] j O^lfTl ). 

N,z 

Let us recall that {Ci}i^j denotes a tiling of by a family of hypercubes of sidelength N. 
We let Jn,z be the set of indices 

Jn,z ■— {j € I, iCj^z Li C(m+l)N,z\C'mN) 7^ 0}- 

and it is clear that the cardinal of Jn^z is We may then write the interaction of 

C(m+i)N,z\CmN with itself as 


Int \C(^m+l)N,z\CmN 1 C*(m+l)JV,^\C*mAf]] 

= - ^ _ g{x-y){dC{x) - dx){dC{y) - dy), 

^ C^Jl ,z*~*^{m+l)iV,2\^miV ) X 32 {m+ 1 ) N ,z \^mN ) 

“I" ^ ^ I'llf [iio',2 L C*(m+l)Af,2\CmAf)) Cj,2 L C*(^_l_l)7V^2\CmA^■]■ 

The previous identity is nothing but writing the interaction of a collection of (possibly truncated) 
hypercubes with itself as the sum of hypercubes-hypercubes interactions plus the sum of self¬ 
interactions. 

Since there are 0(m'^“^) elements in Jn,z we have 


Ep(2) 

N,z 


^{m+l)N,z\^imN\ Cj,z L C*(^m+l)JV,2\^mAf] 


j&JN,z 




It remains to bound the sum of interactions between two disjoint hypercubes. We have, if the 
two hypercubes are not adjacent (since there are only 0(m'^“^) pairs of adjacent hypercubes in 
the sum, all giving a contribution of order 0(1), we may neglect these terms), 


_ g{x - y){dC{x) - dx){dC{y) - dy) 

(Lj'i ,z*~'L^m+l)iV,z \LmW ) X (Lj2 :Z^G ^ mN ) 

> -Cg{disi{Cp,Cj^)), 

with a constant C depending on N and d,s. For any hxed j G J^^z we have (the sum is 
implicitely restricted to non-adjacent hypercubes) 

/ mN 

r~^r^~‘^dr < + 1), 

r^3^JN,z 

with a constant C depending on N and on d,s. The element of volume is only r'^~^ because 
we are summing terms on the boundary of CmN- We may then estimate the sum of pairwise 
hypercube interactions as 

E 5(dist(0,-i,0,2)) < 0(m'^"^)(0(m'^"^"^) 0(1)) = o(m'^), 

jl¥^j2&JN,z 

because we are restricted to the non-Coulomb cases for which s > d — 2. It proves (I5.23p hence 
also (|5.22l) . Combining (|5.22l) and (I5.2ip we obtain (|5.19l) . 
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Study of Inte[CmN\C(rn-i)N,z-:C(rn+i)N,z\CmN]’ It remains to control the last term in 
the right-hand side of (I5.17p . which is another boundary-boundary interaction and is bounded 
as in (15. 2, ip . 

Conclusion for the energy. Finally, combining (|5.16l) and (I5.18p we obtain that 
^p(2) [ln-t[C(^_|_]^)7v 2 ) IrLt[C*mAf) C*mAf]] ^ oijTl ), 

^N,z 

with a o{m'^) depending on N and d,s, but uniform in z G Cj\f. Combining this estimate with 
p5.14p we thus conclude that 

limsuplimsupE (3) [lnt[CmN,CmN]] < 

N—i'OO m^oo N 

Step 3. Entropy and conclusion. The convergence of the entropy follows easily from (15.1211 
and the definition of the relative specific entropy, indeed we have 

ent[pj^^|n] = ^Ent[Piv|nc7^] + o(l). 

(o\ 

To summarize, we have shown that the sequence {P^ }n satisfies the requirements of Propo¬ 
sition 15.31 which concludes the proof. □ 


6 High-temperature limit 

In this section we apply the results of Section [5] to study the limit as /3 ^ 0 (the high-temperature 
limit) of the minimizers of Pg. We prove their convergence to the law of the Poisson point process 
in all cases (HTH), ([O]), ([L3|) as stated in Theorem [2l 


6.1 Specific Pinsker inequality 


The well-known Pinsker inequality gives an upper-bound on the total variation distance between 
probability measures in terms of their Kullback-Leibler divergence: 

(6.1) |P-Q|TV<y^Ent(P||Q) 

where \P — Q\tv is the total variation defined by 

(6.2) \P — QItv := sup{P(74) — Q{A),A measurable}. 


Combining the Pinsker inequality with the property (|2.13l) of the specific relative entropy we 
get for any stationary random point process P the following specific (infinite-volume) Pinsker 
inequality: 


(6.3) sup n^.iTv ^ nt[P|n]. 

N>1 \Cn\^ ^ ^ 

Since the total variation convergence implies weak convergence of probability measures, it 
is clear that the convergence in specific relative entropy sense implies the weak convergence of 
random point processes i.e. if a sequence of stationary random point processes {Pk}k satisfies 


lim ent[Pfc|n] = 0 

fc—>oo 

then the sequence {Pk}k converges to 11. 
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6.2 Finite energy approximation the Poisson point process 

Since the two-point correlation function of the Poisson point process satisfies /92,n = 1 we 
clearly have W™*(II) = 0 in all cases (HH), dLlI) and (11.31) . In the case (11.31) it thus follows 
that is finite, according to the electric-intrinsic inequality of Section 15.11 For the 

one-dimensional Log-gas it has been proven in |LS15] that W®*®'^(n) = -|-oo, and the answer is 
unknown in the two-dimensional Log-gas case (the result of Section 15.11 is not enough because 
piog(n) is infinite). However we may always construct random point processes which converge 
in entropy sense to 11 and whose renormalized energies are finite. 

Lemma 6.1. There exists a sequence {'Kk^k of stationary random point processes in Vs,i{X) 
satisfying 

1. is finite for all k. 

2. limfc^ooent[7rfe|n] = 0. 

Proof. For any A: > 1, let be a tiling of by a countable family of disjoint copies of 

the hypercube Ck, and let {B^}j be the law of a family of independent Bernoulli point processes 
with k'^ points in C^. We let be the random point process consisting of the union of all 
for i G I. Finally we define by averaging over a “fundamental domain” i.e. we let 

(6.4) TTfc := Uk + vTfc, 

where Uk is a uniform random variable in Ck (if C is a point configuration we let x + C denote 
the point configuration {x + p,p G C}, cf. also (ll.lOp i. 

The random point processes dehned this way are clearly stationary and of intensity 1. 
The two-point correlation function of is easy to compute: 

, . . , , I 1 — 7 ^ if X and y belong to the same hypercube Cf, 

(6.5) P2ix,y) = < . 

1 otherwise. 


The two-point correlation function of vr^. could be deduced from (16.5p by averaging p 2 {x.,y) over 
translations of both coordinates by a vector in Ck- Let us simply observe that p 2 {x,y) — 1 is 
bounded (because p\ is) and has compact support (e.g. p2{x,y) = 1 as soon as \x — y\ > Vdk) 
which implies (using the expression (I3.12p i that W“*(7rfc) is finite. Moreover, observing that 


I 

Jcl 


{pi{x,y) - 1 ) = -1 


we also get that 


f 

Jcl 


{P2{x,y)-1) = -1 + 0{R^ ^) 


which implies in view of (15.ip that hence T>^°^{7rk) is zero. Using the 

electric-intrinsic inequality we conclude that W®*®'^(7rfc) is hnite for all A; > 1 which proves the 
first point. 

We are left to prove the second point of the lemma. Let be the measurable subset of point 
configurations which have exactly k^ points in each hypercube C^, and let £k be the subset 
obtained by averaging Ej, over Ck, more precisely we let 


( 6 . 6 ) 


£k 1^ x C 
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where the sum x + C is dehned as above. By definition tt/j coincides with the law of the Poisson 
point process conditioned to the event £k- For any i? > 0 we may thus estimate the relative 
entropy 


(6.7) Ent [(7rfc)c«|n|c^] = - logn|c^(ffe). 

Since 8^^ C £k we may bound below 'n.\^Q^{£k) by which is easier to compute, indeed 

we only have to estimate the probability that 




'rY ^ R'^ 

k 


disjoint hypercubes of sidelength k receive exactly k'^ points, and we can bound below 
by 


n|CH( 4 ) > 



NR,k 


An elementary estimate using Stirling’s formula shows that 

( 6 . 8 ) 

with a universal constant C. We deduce from (16.7|) and (16.8p that 


jjjEnt [(>rt)c,|n|cj = 0(1-“) 

hence by definition of ent[jn] we also have ent[7rfc|II] = 0{k~'^), which proves the second point 
of the lemma. □ 


6.3 Proof of Theorem [2] 

From Lemma l6.II the proof of Theorem [2] is straightforward. 

Proof. For any /I > 0, let P /3 be a minimizer of J-p. In particular we have 

+ ent[P^|n] < + ent^jn] 

for any A: > 1, where {7rk}k is the sequence of random point processes constructed in Lemma 
16.11 Since is bounded below by some constant depending only on d, s we have 

(6.9) sup ent[P/ 3 |n] < ent[7rfc|n] + /3 (w®'®®(7rfc) - minW®'®®) . 

J^I3iPi3)=kainTp ^ 7 

Since ent[7rfc|n] = Ofc(l) and W®'®®(7rfc) is always finite, we get (II.9|) by considering k large enough 
and fd small enough (depending on k). The fact that convergence in entropy sense implies weak 
convergence was observed in Section 16.11 □ 

Since Sine^ was proven to be a minimizer of J-p for the one-dimensional Log-gas, we get 
Corollary 11.21 as an immediate consequence of Theorem [2l This convergence result was recently 
established in |AD14] by analysing the family of coupled diffusion processes dehning the point 
processes Sine^j. Here we rely only on the fact that the Sine^j process minimizes the free energy 
functional Tp. 
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7 Low temperature limit in one dimension 

In this section we prove Theorem [3] i.e. we use the link between and to give a 

minimization result on the energy in the one-dimensional case. As can be expected the minimizer 
of W®*®® is attained by a “crystalline state” which in dimension 1 corresponds simply to the lattice 
Z. In the remaining of this section we deal with the cases dEU) or with d = 1 and 0 < s < 1. 

In |SS14] (in the one-dimensional logarithmic case) the minimality of W(Z) among the ener¬ 
gies of periodic configurations was proven using an explicit formula valid in the periodic setting, 
together with the convexity of the interaction kernel. An argument of approximation by periodic 
configurations was then used to prove that Z is a global minimizer of the energy (however, it is 
not unique). In [Lebl5| we turned this convexity argument into a quantitative estimate in order 
to bound below the difference W®^®®(P) — W®^®®(Pz) in terms of the two-point correlation function 
of P, first in the periodic case, then in the general stationary case using the same kind of ap¬ 
proximation. It was enough to prove that Pz is the unique minimizer of W®*®® among stationary 
point processes in the case d = l,s = 0. It also yields the fact that if W®^®®(Pn) —^ W®^®®(Pz) 
then the two-point correlation function of Pn converges to that of Pz- In the following we 
prove the same result in all one-dimensional cases, first at the level of hyperuniform random 
point processes (which include periodic point processes) then in the general stationary case us¬ 
ing the approximation argument of Proposition 15.31 We also observe that convergence of the 
two-point correlation functions to that of Pz in fact implies weak convergence of the random 
point processes. 

7.1 The k-th neighbor correlation functions 

In the one-dimensional case the two-point correlation function of a stationary random point 
process admits a decomposition as the sum of the A:-th neighbor correlation functions. 

Let P be in Vs^i[X) such that the two-point correlation p 2 exists as a Radon measure in 
M X M. For any k > 1 we define the /c-th neighbor correlation function p 2 ^k by duality, letting 
for any £ C'c(M X M) 


(7.1) 


J ^P2,k 



<p{x,y) + p{y,x) 

x^y£C,y k-th neighbor of x 


In (17.ip if x,y belong to a point configuration C we say that y is the /c-th neighbor of x if x < y 
and C{[x,y]) = k + 1. We will abbreviate “/c-th neighbor of’ by k.n.o. in the formulas. Since P 
is stationary we may see p 2 ,fc as a measure on M by letting P 2 ,k{x) '■= /92,fc(0,x) (in the rest of 
this section we will use the same notation for both interpretations of P 2 ,k)- 

Lemma 7.1. For any k > 1, for any compactly supported, measurable, even function (/? : M i—)• M, 
we have 


(7.2) 


Ep 


x^y£.Cr\[—L/2^L/2]^y k.n.o. x 



^{x)p2,k{x){l 


x/L). 


Proof We use the definition (dD together with a change of variable (x, y) ^ (x — y,x + y) as 
in the re-writing of W™* in Section 13.21 □ 
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7.2 Minimization: the hyperuniform case 

Lemma 7.2. The unique minimizer 0/ among random point processes in 'Phyp(<^) is the 
random point process Pz defined in (11.101) . Moreover for any such P we have 


(7.3) 


~l~oo n-\-oo 




int 


{Pi)>cY, / 


mm 


(x — k)" 

]^s+2 


, 1 P2,k{x), 


with a constant c depending only on s. 


Proof. Since there are a.s. r points in any interval of length r the A:-th neighbor correlation 
function of P is supported in [k — r,k + r]. We may thus write, for any R > 0 

I rR 2 /-R+f- \ 

R J “ \x\)dx = 9{x) ( ^ P2,k - ij (i? - x)dx. 


Let fiR : X ^g{x){R — x). By definition of W™* (see (I3.12p l we get 


(7.4) 




int 


rR /R+r \ 

(P) = hminf / 'ifnix) I ^ p 2 ,k -Ij dx, 

\fc=i / 


and for Pz it is easy to see that the liminf is actually a lim and we have 

rR 


jint 


(7.5) 

Substracting m in dzaD we get 


= lim 
R—^oo 


rU / 

J i’R{x) E — 1 I dx. 


W™*(P) - W™*(Pz) = liminf [ ifnix) [ V p2,fe - ^ ' 

We may re-write the previous expression as 

^+00 /1^“’'! \ 

( 7 . 6 ) W“'(P) - W“*(Pz) = liminf / finix) V {p2,k - 4 ) + 

Jo \ / 


IRi 

-y.^k 


k=l 


where the error term pR^r is bounded using the fact that p 2 ,fc is supported on [k — r,k + r] and 
that IfiRl is decreasing on [R — 2r,R] (for R large enough). 


rR / R+r 

(7.7) \ER^r\ < / f’Rix) y 

■^0 \k=lR- 


[R\ 

P2,k+ 'y 5k\ <Cr\il){R-2r)\= or{1). 

= \H—r\ k=[R—r\ 


Let us now observe that fiR is a convex function, more precisely for x G (0, -|-oo) we have 
f^Rix) > for some positive constant c depending on s. Moreover for all fe > 1, since P is of 
intensity 1 we have f p 2 ,fc = 1 iP 2 ,k is the probability law of the fc-th neighbor) and since P is 
periodic the expectation f xp 2 ,k(x) is finite and thus equal to k (the A:-th neighbor is in average 
at distance k). Combining this observation with the convexity estimate we may write 


r-hoo P+00 

/ 'lpR{x){p 2 ,k - 4 ) > / 

Jo Jo 


(x — kfi 


Q max(x, kfi^^ 


r+co 

P2,k{x) >c min 
Jo 


{x — kfi 

ks+2 


, 1 P2,ki^)^ 
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with c depending only on s. Inserting this bound in (|7.6p and using (17.7p we get 

^ + 00 /(rf._U\2 \ 

Wmt(p) _ p J ^ ^ 1 j P2,k{^) + 0^(1) 

Hence finally by taking the limit ii —)• oo we obtain that 

(7.8) w“*(p) - w“*(Pz) > -°2,fc(x), 

which proves (US]). This lower bound implies that Pz is a minimizer of among hyperuniform 
random point processes and also that it is unique. Indeed W™*(P) = W™*(Pz) implies (by (I7.3p i 
that p 2 ,k = for all A; > 1, hence the two-point correlation function of P coincides with the one 
of Pz, which is enough to conclude that P = Pz (see e.g. |Lebl5] l. □ 

7.3 Proof of Theorem [3] 

We may now give the proof of our minimization result for First let us observe that 

W®i®c(Pz) = W™*(Pz). The inequality < is true by Proposition 15.11 Moreover Proposition 
15.31 ensures that there exists a sequence of hyperuniform random point processes such that 
limsupjv^oo W‘“*(PAr) < W®^®®(Pz). By Lemma [72] above we know that W™^{Pn) > W“*(Pz), 
hence in fact W‘“*(Pz) < W®^®®(Pz) and equality holds. 

Proof. Step 1. Minimization o/ W®*®®. Let P be a minimizer of W®*®® on Vs^i{X). From 
Proposition 15.31 we get a sequence {Pn}n of hyperuniform random point processes converging 
to P and such that {W™*(Pn)}„ converges to W®^®®(P). For any k,n> I let denote the 
k-th. neighbour correlation function of Pn- In the hyperuniform case, (|7.3I) implies that for any 
M > 0 and A: > 1 we have 

(7.9) W-*(P„) - W“*(Pz) > 

and the right-hand side is bounded below by (see (I7.2p i 

(pk{x-y) 

x,yGCn[—M^M],y k.n.o. x 

where we let <Pkix) ■= min , 1^- Since {Pn\n converges to P we have 


r2M 

/ min 

Jo 


{x — ky^ 

ks+2 




X \ 

~ ~ 


= 


(7.10) Ep„ 

1 

1 

9- 

W 

1_ 

= Ep 

1 

1 

..se 

9- 

W 

1_ 


Cn[—M,M],y k.n.o. x 


Cn[—M,M],y k.n.o. x 


r2M 

/ ^kix)p2,ki^) 

Jo 


1 - 


Combining (|7.9p and ()7.10p we see that 


rM 

/ min 

Jo 


(x — ky 

ks+2 


2M 


, 1 ) P 2 ,fc(x) < 2 (W“*(P„) - W^"'(Pz)) + On(l), 


+ On(l) 

1 /■“ 

j + On(l) > 2 J ^kP2,k{x) + On(l). 
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but as n —>• oo we have —)• = W™*(Pz)- It implies that 


^ ,1^ P2,fc(3;) = 0 

for all M > 0 and fc > 1. Finally we get that p 2 ,A: = for all A: > 1 and we conclude as in the 
proof of Lemma 17.21 that P = Pz, which ensures that Pz is the unique minimizer of on 

Ps,lW- 

If P is not a minimizer the same argument shows that 



as in the hyperuniform setting. 

Step 2. Energy of minimizers of Pp tends to On the other hand we claim that 

if {F’/ 3 }/ 3 >o is a family of minimizers of Pfs then we must have 

(7.12) lim 

/3^.0 

To prove (I7.12D we cannot directly evaluate over Pz and compare with Pp because ent[Pz|n] 
is infinite. However we may argue as in Section 16.21 and show that there exists a sequence {vr^lfc 
of stationary random point processes in Vs^i[X) satisfying 

1. ent[7rfc|n] is finite for all k. 

2. hmfc^ooW®'®®[7rfc] =W®i®®(Pz). 

Such a sequence can be constructed by chosing a “vibrating” approximation of Pz- For any 
k > 1 we let {Vk^rn}mez be a countable family of i.i.d. random variables distributed uniformly 
in then we let vr^ be the random point process 


'^k ^rn+Vk,m 

meZ 


and hnally we define vr^ by averaging vr^ over [0,1]. It is easy to check that vTfc is a stationary 
random point process of intensity 1. In fact may equivalently be defined as a renewal process 
with increments distributed as 1 + 14^2 ~ 14,i- The specific relative entropy of vr^ coincides with 
its “entropy rate”, it is finite (see |DVJ08l Section 14.8.]) and blows up as A: —?■ oo like the entropy 
of 14,2 — 14 , 1 - Concerning the energy, we have the bound 

W®'"®(Pz) < W®^®®(7rfc) < W“*(7rfe) 


and the fact that W™*(7rfc) converges to W™*(Pz) = W®*®®(Pz) can be checked directly with 
the help of the formula defining W“*. Indeed the two-point correlation function of vr^ may be 
written as 

P2,TTk = ^fJk{m + -) 
m^'L 

where is a triangular “hat function” of width ^ and integral 1. For any m G Z and 7? > 0, a 
mean value argument shows that 



m + ■) log |x|(l 


\x\ 

"r 


) — log |m|(l 


R ~ kf mR j 
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Consequently, we get 


lim snp 

R^oo 



P2,z) (1 




and we obtain that lim^^oo W™*(7rfc) = 

With the help of the sequence {vr^j^ we obtain (I7.12|) by arguing as in Section lUlBl indeed 
since ent[P^|II] is always non-negative we have 

+ ^ent[7rfc|n] < + ^ent^ln] + 0^(1) 

and ()7.12p follows by chosing k, /3 large enough. 

Step 3. Convergence of the two-point function of minimizers of Pg. Combining ()7.1ip and 
(I7.12P we see that if denotes the k-ih neighbor correlation function of Pg, we have for any 


X 



min 


/ (x — k)"^ 
\ A :®+2 


10 


as /3 ^ oo. Arguing as in the proof of |Lebl5[ Lemma 2.3] we deduce that converges to 
in the distributional sense as /3 —>• oo. Let us observe that 


P2,z := ^ hk 
kez* 

is the two-point correlation function of Pz- 

Step 4. Weak convergence of the minimizers ofpg. It is not hard to see that this convergence 
implies in fact the weak convergence of Pg to Pz as ft ^ oo. For any ^ > e > 0 let be a 
smooth non-negative function which is equal to 1 on the set Uk£z[k — t + £,k — e] and to 0 
on Z. For any T > 0 we let be a non-negative continuous function such that tt^s = 1 on 
[—T + £,T — e] and 0 outside [—T, Tj. We let ipe^T be the continuous, compactly supported map 


P‘e,T{x,y) := Xeix - y)TT,e{x)TT,eiy) ■ 

Let Ax^e be the event “there is no pair (x, y) of points of the configuration in [—T P e,T — e] 
such that |x — y| € Uk^z[k — 1 + e,k — ej”. Since f Pe,TP 2 ,'L = 0 and since converges to p 2 ,z 
as /3 —>■ oo we have 


(7.13) 

as ft ^ oo. In other words, with probability tending to 1 as /3 —>■ oo, a configuration under 
Pp locally looks like a (translated) subset of Z in which all the points have been displaced at a 
distance at most e. 

The variance under Pp of the number of points in [—T, T] is bounded as /3 —)• oo, because it 
is controlled by the energy W®^®‘^(P^), which itself converges (to W®*'^'^(Pz)). This follows from 
the discrepancy estimates (see e.g. |Lebl5[ Lemma 2.1] or [LSISI Lemma 3.10]). In particular we 
have uniform integrability under Pp of the number of points in [—T, T] as /3 —)• oo. In particular, 
conditioning Pp to we have an average of 2T -|- o(l) points in [—T, T] as ft ^ oo, because 
(TmH holds. 

Finally we deduce that for any e > 0, with probability tending to 1 as /3 ^ oo, a conhguration 
under Pp locally looks like a translate of Z in which all the points have been displaced at a 
distance at most e. This implies the convergence of Pp to Pz as ft ^ oo. 

□ 
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